PROCEEDINGS 

°F SCIENCE 



Phase diagram of QCD with two degenerate 
staggered quarks 



o 

(N 

o 
O 

00 



> 
o 

cn 

d 



p. Cea 

Dipartimento di Fisica dell'Universita di Ban, 1-70126 Ban, Italy and INFN, Sezione di Bari, 
1-70126 Bari, Italy 



E-mail: paolo . ceaSba . inf n . it 



L. Cosmai* 

INFN, Sezione di Bari, 1-70126 Bari, Italy 
E-mail: 



leonardo . cosmaiSba . inf n . it 



M. D'Elia 

Dipartimento di Fisica dell'Universita di Genova, 1-16146 Genova, Italy and INFN, Sezione di 

Genova, 1-16146 Genova, Italy 

E-mail: jnassimo . delia@ge . inf n . it 



A. Papa 

Dipartimento di Fisica dell ' Universitd della Calabria, 1-87036 Rende ( Cosenza ), Italy and 
INFN, Gruppo collegato di Cosenza, 1-87036 Rende ( Cosenza), Italy 



E-mail: papaScs . inf n . it 



F. Sanfilippo 

Dipartimento di Fisica, Universitd di Roma "La Sapienza " and INFN, Sezione di Roma, 
Piazzale A. Moro 5, 00185 Roma, Italy 
E-mail: 



francesco. sanf ilippoSromal . inf n . it 



We present preliminary results about the critical line of QCD with two degenerate staggered 
quarks at nonzero temperature and chemical potential, obtained by the method of analytic contin- 
uation. As in our previous studies with different numbers of colors and flavors, we find deviations 
from a simple quadratic dependence on the chemical potential. We comment on the shape of the 
critical line at real chemical potential and give an estimate of the curvature of the critical line, 
both for quark chemical potential and isospin chemical potential. 
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1. Introduction 

Many important issues, such as heavy ion collisions, evolution of the early universe, physics 
of the compact stars, are related to the knowledge of the QCD phase diagram. Understanding the 
physical phases of QCD [1] in ordinary as well as in extreme environment of high temperature and 
high baryon number is therefore a major research goal. Lattice QCD is the main nonperturbative 
tool to investigate the QCD phase diagram. Unfortunately lattice QCD simulations at finite quark 
density are plagued by the well-known sign problem. As a matter of fact, importance sampling 
requires positive weights in the partition function 



but when H ^0 this is no longer true, indeed det[M(jU)] becomes complex in the case of SU(3) 
theory and finite quark chemical potential. There are few exceptions where det[M(/i)] is real (pos- 
itive for an even number of flavors), namely for two-color QCD, for isospin chemical potential 
and for imaginary values of the quark chemical potential. The sign problem can be addressed by 
means of several techniques, each one, however, suffering from limitations. In the last few years 
we approached the sign problem by the method of analytic continuation [2^]. To this purpose we 
started our investigations [5-7] from some special cases where Monte Carlo simulations are both 
feasible at real and at imaginary chemical potential. The reason was to have the possibility to test 
the results of the analytical continuation by means of direct simulations at real chemical potential. 
In the following sections we review some of our earlier results and give preliminary data for the 
critical line of QCD with Nf = 2 standard staggered fermions. 

2. Investigations about analytic continuation 

The first lesson we learned in studying the phase diagram of two-color QCD [6] was that non- 
linear terms in the dependence of jSc on in general cannot be neglected. Indeed the prediction 
for the pseudocritical couplings at real chemical potentials may be wrong if data at imaginary 
are fitted according to a linear dependence in n^. In Fig. l(Left) the discrepancy between the linear 
extrapolated critical line and the direct determinations of ^c{i^^) is shown for the case of two-color 
QCD with Nf = 8 degenerate staggered fermions. In the case of SU(3) isospin chemical potential 
with Nf = % degenerate staggered fermions [7], instead, few interpolations can be found which 
correctly detect deviations from the linear behavior in at imaginary chemical potentials and 
lead to consistent extrapolations to real /i, in agreement with the direct determinations there. In 
Fig. I (Right) we give an example of a good fit (and the related extrapolation) of the critical line, 
using a ratio of polynomials ^{}x'^) / i)'{}X^). 

We have also revisited [8] the determination of the pseudocritical line of QCD with Nf = 4 
degenerate staggered quarks at nonzero temperature and quark density by the method of analytic 
continuation. The sign problem in this case prevents us from performing simulations at real quark 
chemical potential. Notice that in this case the linear fit is not able reproduce the data for pseudo- 
critical couplings even at imaginary quark chemical potentials. To derive the critical line we ex- 
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ploited analytical continuation by means of polynomial or ratio of polynomials fit: 

In Fig. 2(Left) we display the interpolation of the critical line at imaginary values of the quark 
chemical potential obtained by mean of a ^{iJ.^) / ratio of polynomials fit. 

A nice alternative to the polynomial fits Eq. ( p] ) is given by what we call "physical fit". The 
idea is to write the interpolating function in physical units and to deduce from it the functional 
dependence of jSc on /x^, after establishing a suitable correspondence between physical and lattice 
units. The natural, dimensionless variables of our theory are 7/7^(0), where ^^(O) is the critical 
temperature at zero chemical potential, and n/T. The ratio T /Tc{0) is deduced from the relation 
T = l/{Nta{j5)), where Nt is the number of lattice sites in the temporal direction and a(j8) is the 
lattice spacing at a given j3 (we discard here the dependence on the bare quark mass am) . We 
use for fl(j8) the perturbative 2-loop expression with the given number of colors and flavors. Now, 
adopting the 3 -parameter function 



r,(o) 



H 2 



l+CAtVr,'(M) 



l+AAtV7:2(At)+BAtV7;4(At) 



(2.2) 



we are led to the following implicit relation between jSc and /x^ 



In Fig. 2(Right) we display the interpolation of the critical line at imaginary values of the quark 



chemical potential obtains by mean of the 2-loop "physical fit" of Eq. ( |2.3[ ). 

Regrettably, even if we are able to have successful fits of the data at imaginary chemical 
potential, the results of the extrapolation of these fits to real values of the quark chemical potential 
are quite consistent only up to /i/r ~ 0.6 (see Fig. ^, where our results are compared with data in 
the literature collected in Ref. [9]). 
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Figure 2: SU(3), Nf = 4, finite quark chemical potential. Fits to the critical couplings: ratio of a 4th- to 
2th-order polynomial (Left) and 2-loop "physical fit" according to Eq. (2.3). 
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Figure 3: Comparison of our extrapolations with other determinations in the Uterature. For the sake of 
readabihty, our extrapolations have been plotted without error bands. Legenda: D'Elia, Lombardo, Ref. [4]; 
Azcoiti et al., Ref. [10]; Fodor, Katz, Ref. [11]; Kratochvila, de Forcrand, Ref. [12]. 



3. QCD with two degenerate staggered quarks: preliminary result 

In this Section we present preliminary results for the critical line of QCD with two degenerate 
staggered flavors. We are going to consider the two cases of finite quark chemical potential and of 
finite isospin chemical potential. 

We performed simulations using standard action with staggered fermion mass am = 0.05 on 
a 16^ X 4 lattice. Simulations were done using the RHMC algorithm, properly modified for the 
inclusion of a finite chemical potential. We collected typical statistics of 10 thousand trajectories. 
The critical coupling at given chemical potential, ^dp-^), is determined as the value for which the 
susceptibility of the chiral condensate exhibits a peak. This value can be deduced by means of a 
Lorentzian fit to the peak or through a reweighting analysis. 

We learned from our previous studies discussed in Section 2 that non-linear terms in /i^ cannot 
be neglected. We also learned that good fits are achieved by means of two kinds of interpolations. 
The first one is given by polynomial or ratio of polynomials (see Eq. (p])). The second one 
consists in writing down the interpolating function in physical units and in performing a fit to the 



data through this implicit relation between jSc and }l (see Eq. (2.3)). Actually, the linear fit to the 
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critical couplings versus imaginary values of the quark (or isospin) chemical potential has very bad 
quality and we are led to exploit, even for SU(3) and Nf = 2 staggered flavors, the functional forms 
given by Eq. (2J_) or by Eq. (23). In Fig. ^ we give a sample of good interpolations obtained to 
data for quark chemical potential and for isospin chemical potential by means of the ratio of 4-th 



order polynomial to 2-nd order polynomial or employing the "physical fit" Eq. (2.3). 
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Figure 4: (Left) SU(3) Nf — 2, finite quark chemical potential. Fit to the critical line by means of ratio of 
4th order to 2nd order polynomial (top) or by means of the "physical fit" Eq. (2.3) (bottom). (Right) SU(3) 
Nf = 2 isospin chemical potential. Ratio of polynomial fit (top) and "physical" fit (bottom). 



The circumstance that there are significant deviations from a simple linear behavior in is 
evident by looking at Fig. 5(Left), where a comparison between the ratio of polynomials fit and the 
linear fit to data at small imaginary chemical potential are displayed. The linear fit is restricted to 
imaginary chemical potential values such that ;^^/dof < 1. This corresponds to imaginary chemical 
potential values with /i^ < 0.1. 

In Fig. ||(Left) we show a comparison between the extrapolations to real values of the quark 
chemical potential starting respectively from three of the more successful interpolations of imag- 
inary quark chemical potential data. If we consider the extrapolation done starting from the 
"physical" fit to the quark chemical potential data, and in the hypothesis that we can trust it 
down to the T = axis we can give an estimate of the critical value of at 7 = 0. We find 
jU(r = 0) = 3.284(64)rc(0) to be compared with ix{T = 0) = 2.73(58)7^(0) obtained in Ref. [13] 
with Nf = 2 Wilson fermions. 

However, while two of the extrapolations are quite consistent with each other, the third one 
(the sixth order constrained, obtained with a sixth order polynomial fit where the first term is 
constrained to the value obtained by a linear fit to small imaginary chemical potential) differs from 
the other two for values iJ./{nT) larger than approximately 0.2. In Fig. ^Right) we compare the 
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Figure 5: (Left) SU(3) Nf — 2, comparison between the ratio of polynomials fit (red lines) and the linear 
fit to data at small imaginary chemical potential (blue lines). (Right) Linear extrapolations starting from 
small values of the imaginary quark chemical potential (red lines) and small values of the imaginary isospin 
chemical potential (blue lines). The full dots are the results of direct numerical simulations at real values of 
the isospin chemical potential. 
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Figure 6: (Left) SU(3) Nf ^2 isospin chemical potential. Comparison between ratio of 4th order to 2nd 
order polynomial fit and linear fit to small values of the imaginary isospin chemical potential. (Right) Ex- 
trapolations from linear fits at small chemical potential both for the case of isospin chemical potential (blue 
lines) and quark chemical potential (red lines). The full dots refer to data obtained from direct simulations 
at real isospin chemical potential. 



same extrapolations in the case of isospin chemical potential. The outcome is almost the same. The 
full circles are the values of the critical coupling from direct simulations at real isospin. The first 
two data are in agreement with all three different kinds of extrapolation. The third point here is in 
between. We are now running simulations at larger values of real isospin chemical potential. 

Finally we are able to give an estimate of the curvature of the critical line at = defined as 
the derivative d^^ [pL^) / dpL^ evaluated at = 0. By means of the extrapolations from the linear fits 
at small chemical potential (see Fig. ^Right)) 



1 + 



«i i-^) (3.1) 



Tc{Q) \kTJ 
we get the following values for the isospin and for the quark chemical potential: 

fli = —0.470(13) isospin chemical potential (3.2) 

a\ =—0.522(10) quark chemical potential . (3.3) 
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The curvatures in the two cases differ for about 5 standard deviations. The curvature of the critical 
line for quark chemical potential is in good agreement with ai = —0.500(34) from Ref. [14], where 
am = 0.025. 



4. Conclusions 

We have applied the method of analytic continuation to study the pseudo-critical line jS^ (/i^) in 
QCD with Nf = 2 degenerate staggered fermions in the case of non-zero quark chemical potential 
and in the case of non-zero isospin chemical potential. We have estimated the curvature of the 
critical line. As in our previous investigations we find significant deviations from a simple linear 
behavior in /i^. Notwithstanding, there are several kinds of functions able to interpolate the critical 
line at imaginary values of the chemical potential, leading to extrapolations which diverge from 
each other at large real jx. 
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